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Abstract. In this paper we present and discuss for the first time a large data set of UBVRI night sky brightness 
measurements collected at ESO-Paranal from April 2000 to September 2001. A total of about 3900 images obtained 
on 174 different nights with FORS1 were analysed using an automatic algorithm specifically designed for this 
purpose. This led to the construction of an unprecedented database that allowed us to study in detail a number 
of effects such as differential zodiacal light contamination, airmass dependency, daily solar activity and moonlight 
contribution. Particular care was devoted to the investigation of short time scale variations and micro-auroral 
events. The typical dark time night sky brightness values found for Paranal are similar to those reported for other 
astronomical dark sites at a similar solar cycle phase. The zenith-corrected values averaged over the whole period 
are 22.3, 22.6, 21.6 20.9 and 19.7 mag arcsec -2 in U, B, V, R and / respectively. In particular, there is no evidence 
of light pollution either in the broadband photometry or in the high-airmass spectra we have analysed. Finally, 
possible applications for the exposure time calculators are discussed. 

Key words, atmospheric effects - site testing - light pollution - techniques: photometric 



1. Introduction 

The night sky brightness, together with number of clear 
nights, seeing, transparency, photometric stability and 
humidity, are some of the most important parameters 
that qualify a site for front-line ground-based astronomy. 
While there is almost no way to control the other char- 
acteristics of an astronomical site, the sky brightness can 
be kept at its natural level by preventing light pollution 
in the observatory areas. This can be achieved by means 
of extensive monitoring programmes aimed at detecting 
any possible effects of human activity on the measured 
sky brightness. 

For this purpose, we have started an automatic survey 
of the UBVRI night sky brightness at Paranal with the 
aim of both getting for the first time values for this site and 
building a large database. The latter is a fundamental step 
for the long term trend which, given the possible growth 
of human activities around the observatory, will allow us 
to check the health of Paranal's sky in the years to come. 

The ESO-Paranal Observatory is located on the top of 
Cerro Paranal in the Atacama Desert in the northern part 
of Chile, one of the driest areas on Earth. Cerro Paranal 
(2635 m, 24° 40' S, 70° 25' W) is at about 108 km S of 
Antofagasta (225,000 inhabitants; azimuth 0?2), 280 km 
SW from Calama (121,000 inhabitants; azimuth 32?3), 152 
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km WSW from La Escondida (azimuth 32?9), 23 km NNW 
from a small mining plant (Yumbes, azimuth 157.°7) and 
12 km inland from the Pacific Coast. This ensures that the 
astronomical observations to be carried out there are not 
disturbed by adverse human activities like dust and light 
from cities and roads. Nevertheless, a systematic monitor- 
ing of the sky conditions is mandatory in order to preserve 
the high site quality and to take appropriate action, if the 
conditions are proven to deteriorate. Besides this, it will 
also set the stage for the study of natural sky brightness 
oscillations, both on short and long time scales, such as 
micro-auroral activity, seasonal and sunspot cycle effects. 

The night sky radiation has been studied by several au- 
thors, starting with the pioneering work by Lord Rayleigh 
in the 1920s. For thorough reviews on this subject the 
reader is referred to the classical textbook by Roach & 
Gordon ( 1973) and the recent extensive work by Leinert 
et al. l|1998fl , which explore a large number of aspects con- 
nected with the study of the night sky emission. In the 
following, we will give a short introduction to the subject, 
concentrating on the optical wavelengths only. 

The night sky light as seen from ground is generated by 
several sources, some of which are of extra-terrestrial na- 
ture (e.g. unresolved stars/galaxies, diffuse galactic back- 
ground, zodiacal light) and others are due to atmospheric 
phenomena (airglow and auroral activity in the upper 
Earth's atmosphere). In addition to these natural compo- 
nents, human activity has added an extra source, namely 
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Fig. 1. Night sky spectrum obtained at Paranal on 
February 25, 2001 02:38UT in the spectral region covered 
by B, V, R and / passbands (from top to bottom). The 
original FORS1 1800 seconds frame was taken at 1.42 air- 
masses with a long slit of 1" and grism 1501, which pro- 
vide a resolution of about 22 A (FWHM). The dashed 
lines indicate the passband response curves. Flux calibra- 
tion was achieved using the spectrophotometric standard 
star Feige 56 (Hamuy et al. 1992) observed during the 
same night. The absence of an order sorting filter probably 
causes some second order overlap at wavelengths redder 
than 6600 A. 

the artificial light scattered by the troposphere, mostly in 
the form of Hg-Na emission lines in the blue-visible part 
of the optical spectrum (vapour lamps) and a weak con- 
tinuum (incandescent lamps). While the extra-terrestrial 
components vary only with the position on the sky and 
are therefore predictable, the terrestrial ones are known 
to depend on a large number of parameters (season, geo- 
graphical position, solar cycle and so on) which interact 
in a largely unpredictable way. In fact, airglow contributes 
with a significant fraction to the optical global night sky 
emission and hence its variations have a strong effect on 
the overall brightness. 

To illustrate the various processes which contribute to 
the airglow at different wavelengths, in Figure ^ we have 
plotted a high signal-to-noise, flux calibrated night sky 
spectrum obtained at Paranal on a moonless night (2001, 
Feb 25) at a zenith distance of 45° , about two hours after 
the end of evening astronomical twilight. In the B band 
the spectrum is rather featureless and it is characterised 
by the so called airglow pseudo-continuum, which arises in 
layers at a height of about 90-100 km (mesopause). This 
actually extends from 4000A to 7000A and its intensity is 
of the order of 3xl0~ 7 erg s _1 cm~ 2 A -1 sr _1 at 4500A. 



All visible emission features, which become particularly 
marked below 4000A and largely dominate the U pass- 
band (not included in the plots), are due to Herzberg and 
Chamberlain O2 bands (Broadfoot & Kendall 1968). In 
light polluted sites, this spectral region is characterised 
by the presence of Hg I (3650, 3663, 4047, 4078, 4358 and 
5461A) and Nal (4978, 4983, 5149 and 5153A) lines (see 
for example Osterbrock & Martel, 1992) which are, if any, 
very weak in the spectrum of Figure ^ (see Sec. El for a 
discussion on light pollution at Paranal). Some of these 
lines are clearly visible in spectra taken, for example, at 
La Palma (Benn & Ellison 1998, Figure 1) and Calar Alto 
(Leinert et al. 1995, Figs. 7 and 8). 

The V passband is chiefly dominated by [OI]5577A 
and to a lesser extent by Nal D and [OI]6300,6364A dou- 
blet. In the spectrum of Figure the relative contribu- 
tion to the total flux of these three lines is 0.17, 0.03 
and 0.02, respectively. Besides the aforementioned pseudo- 
continuum, several OH Meinel vibration-rotation bands 
are also present in this spectral window (Meinel 1950); 
in particular, OH(8-2) is clearly visible on the red wing 
of Nal D lines and OH(5-0), OH(9-3) on the blue wing 
of [OI]6300A. All these features are known to be strongly 
variable and show independent behaviour (see for example 
the discussion in Benn & Ellison 1998), probably due to 
the fact that they are generated in different atmospheric 
layers (Leinert et al. 1998 and references therein). In fact, 
[OI]5577A, which is generally the brightest emission line in 
the optical sky spectrum, arises in layers at an altitude of 
90 km, while [OI]6300,6364A is produced at 250-300 km. 
The OH bands are emitted by a layer at about 85 km, 
while the Na ID is generated at about 92 km, in the so 
called Sodium-layer which is used by laser guide star adap- 
tive optic systems. In particular, [OI]6300,6364A shows a 
marked and complex dependency on geomagnetic latitude 
which turns into different typical line intensities at dif- 
ferent observatories (Roach & Gordon Il973|l . Moreover, 
this doublet undergoes abrupt intensity changes (Barbier 
1957); an example of such an event is reported and dis- 
cussed in Sec. 

In the R passband, besides the contribution of Nal D 
and [OI]6300,6364A, which account for 0.03 and 0.10 of 
the total flux in the spectrum of Figure ^ strong OH 
Meinel bands like OH(7-2), OH(8-3), OH(4-0), OH(9-4) 
and OH(5-l) begin to appear, while the pseudo-continuum 
remains constant at about 3xl0 -7 erg s _1 cm~ 2 A -1 sr _1 . 
Finally, the / passband is dominated by the Meinel bands 
OH(8-3), OH(4-0), OH(9-4), OH(5-l) and OH(6-2); the 
broad feature visible at 8600-8700A, and marginally con- 
tributing to the I flux, is the blend of the R and P branches 
of O 2 (0-l) (Broadfoot & Kendall 1968). 

Several sky brightness surveys have been performed at 
a number of observatories in the world, most of the time 
in B and V passbands using small telescopes coupled to 
photo- multipliers. A comprehensive list of published data 
is given by Benn & Ellison A U authors agree on the 

fact that the dark time sky brightness shows strong vari- 
ations within the same night on the time scales of tens of 
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minutes to hours. This variation is commonly attributed 
to airglow fluctuations. Moreover, as first pointed out by 
Rayleigh ljP32gj) . the intensity of the [OI]5577A line de- 
pends on the solar activity. Similar results were found for 
other emission lines (Nal D and OH) by Rosenberg & 
Zimmerman ()1967fl . Walker H1988bj) found that B and V 
sky brightness is well correlated with the 10.7 cm solar 
radio flux and reported a range of ~ 0.5 mag in B and 
V during a full sunspot cycle. Similar values were found 
by Krisciunas H1990J) . Leinert et al. l|1995p and Mattila et 
al. I|1996[) . so that the effect of solar activity is commonly 
accepted (Leinert et al. 1998). For this reason, when com- 
paring sky brightness measurements, one should also keep 
in mind the time when they were obtained with respect 
to the solar cycle, since the difference can be substantial. 
A matter of long debate has been the so-called Walker ef- 
fect, named after Walker l|1988b|l . who reported a steady 
exponential decrease of ~0.4 mag in the night sky bright- 
ness during the first six hours following the end of twilight. 
This finding has been questioned by several authors. We 
address this issue in detail later (Sec. Eland Appendix IT))l . 

Here we present for the first time U BVRI sky bright- 
ness measurements for Paranal, obtained on 174 nights 
from 2000 April 20 to 2001 September 23 which, to our 
knowledge, makes it the largest homogeneous data set 
available. Being produced by an automatic procedure, this 
data base is continuously growing and it will provide an 
unprecedented chance to investigate both the long term 
evolution of the night sky quality and to study in detail 
the short time scale fluctuations which are still under de- 
bate. 

The paper is organized as follows. After giving some in- 
formation on the basic data reduction procedure in Scc.[5J 
in Sec. 13 we discuss the photometric calibration and error 
estimates, while the general properties of our night sky 
brightness survey are described in Sec. 0] The results ob- 
tained during dark time are then presented in Sec. and 
the short time-scale variations are analysed in Sec. El In 
Sec. we compare our data obtained in bright time with 
the model by Krisciunas & Schaefer l|1991f) for the effects 
of moonlight, while the dependency on solar activity is in- 
vestigated in Sec. [HI In Sec. El we discuss the results and 
summarize our conclusions. Finally, detailed discussions 
about some of the topics are given in Appendices A-D. 

2. Observations and basic data reduction 

The data set discussed in this work has been obtained 
with the FOcal Reducer/low dispersion Spectrograph 
(hereafter FORS1), mounted at the Cassegrain focus of 
ESO-Antu/Melipal 8.2m telescopes (Szeifert 2002). The 
instrument is equipped with a 2048x2048 pixels (px) 
TK2048EB4-1 backside thinned CCD and has two re- 
motely exchangeable collimators, which give a projected 
scale of 0'.'2 and 0"1 per pixel (24/im x 24/im). According 
to the used collimator, the sky area covered by the de- 
tector is 6f8x6.'8 and 3.'4x3f4, respectively. Most of the 
observations discussed in this paper were performed with 



Table 1. Typical background count rates expected in 
FORS1 (SR) images during dark time and at zenith. The 
last column reports the time required to have a back- 
ground photon shot noise three times larger than the max- 
imum RON (6.3 e~), which corresponds to a 5% contri- 
bution of RON to the global noise. 



Passband 


Count Rate 


tz 




(e" px" 1 s- 1 ) 


(s) 


U 


0.5 


714 


B 


3.8 


94 


V 


15.8 


23 


R 


26.7 


13 


I 


32.1 


11 



the lower resolution collimator, since the higher resolu- 
tion is used only to exploit excellent seeing conditions 
(FWHM<0"4). 

In the current operational scheme, FORS1 is offered 
roughly in equal fractions between visitor mode (VM) and 
service mode (SM). While VM data are immediately re- 
leased to the visiting astronomers, the SM data are pro- 
cessed by the FORS-Pipeline and then undergo a series of 
quality control (QC) checks before being delivered to the 
users. In particular, the imaging frames are bias and flat- 
field corrected and the resulting products are analysed in 
order to assess the accuracy of the flat-fielding, the image 
quality and so on. The sky background measurement was 
experimentally introduced in the QC procedures starting 
with April 2000. Since then, each single imaging frame ob- 
tained during SM runs is used to measure the sky bright- 
ness. During the first eighteen months of sky brightness 
monitoring, more than 4500 frames taken with broad and 
narrow band filters have been analysed. 

As already mentioned, all imaging frames are automat- 
ically bias and flat field corrected by the FORS pipeline. 
This is a fundamental step, since in the case of imaging, 
the FORS1 detector is readout using four amplifiers which 
have different gains. The bias and flat field correction re- 
move the four-port structure to within — 1 electron. This 
has to be compared with the RMS read-out noise (RON) , 
which is 5.5 and 6.3 electrons in the high gain and low gain 
modes respectively. Moreover, due to the large collecting 
area of the telescope, FORS1 imaging frames become sky 
background dominated already after less than two min- 
utes. The only significant exception is the U passband, 
where background domination occurs after more than 10 
minutes (see also Tabled). The dark current of FORS1 
detector is -2.2 x 10~ 3 s" 1 px" 1 ( Szeifert E?U0"0"|> and 
hence its contribution to the background can be safely 
neglected. 

Since the flat fielding is performed using twilight sky 
flats, some large scale gradients are randomly introduced 
by the flat fielding process; maximum peak-to-peak resid- 
ual deviations are of the order of 6%. Finally, small scale 
features are very well removed, the only exceptions being 
some non-linear pixels spread across the detector. 
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The next step in the process is the estimate of the sky 
background. Since the science frames produced by FORS1 
are, of course, not necessarily taken in empty fields, the 
background measurement requires a careful treatment. 
For this purpose we have designed a specific algorithm, 
which is presented and discussed in Patat l|2002[l . The 
reader is referred to that paper for a detailed description 
of the problem and the technique we have adopted to solve 
it. 

3. Photometric calibration and global errors 

Once the sky background I sky has been estimated, the 
flux per square arcsecond and per unit time is given by 
fsky = Isky/{tex P P 2 ) , where p is the detector's scale (arc- 
sec pix -1 ) and t exp is the exposure time (in seconds). The 
instrumental sky surface brightness is then defined as 

m sky = -1.086 ln(I sky ) + 2.5 log(p 2 t exp ) (1) 

with m sky expressed in mag arcsec -2 . Neglecting the 
errors on p and t exp , one can compute the error on the sky 
surface brightness as 5b sky — 5i sky /I sky . This means, for 
instance, that an error of 1% on I sky produces an uncer- 
tainty of 0.01 mag arcsec -2 on the final instrumental sur- 
face brightness estimate. While in previous photoelectric 
sky brightness surveys the uncertainty on the diaphragm 
size contributes to the global error in a relevant way (see 
for example Walker 1988b), in our case the pixel scale is 
known with an accuracy of better than 0.03 % (Szeifert 
2002) and the corresponding photometric error can there- 
fore be safely neglected. 

The next step one needs to perform to get the fi- 
nal sky surface brightness is to convert the instrumental 
magnitudes to the standard UBVRI photometric system. 
Following the prescriptions by Pilachowski et al. (^989), 
the sky brightness is calibrated without correcting the 
measured flux by atmospheric extinction, since the effect 
is actually taking place mostly in the atmosphere itself. 
This is of course not true for the contribution coming from 
faint stars, galaxies and the zodiacal light, which however 
account for a minor fraction of the whole effect, airglow 
being the prominent source of night sky emission in dark 
astronomical sites. The reader is referred to Krisciunas 
(1990) for a more detailed discussion of this point; here 
we add only that this practically corresponds to set to 
zero the airmass of the observed sky area in the calibra- 
tion equation. Therefore, if M sky , m sky are the calibrated 
and instrumental sky magnitudes, M*, to* are the corre- 
sponding values for a photometric standard star observed 
at airmass z„, and k is the extinction coefficient, we have 
that M sky = (m sky - to*) + K z* + M* + 7 (C sky - C»). 
This relation can be rewritten in a more general way as 
M sky — too + m sky + 7 C sky , where mo is the photo- 
metric zeropoint in a given passband and 7 is the colour 
term in that passband for the color C sky . For example, 
in the case of B filter, this relation can be written as 
B sky = B + b sky + j§_ v (B — V) sky . 
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Fig. 2. UBVRI photometric zeropoints for FORS1 during 
the time range covered by sky brightness measurements 
presented in this work (vertical dotted lines). The thick 
segments plotted on the lower diagram indicate the pres- 
ence of sky brightness data, while the arrows in the up- 
per part of the figure correspond to some relevant events. 
A: water condensation on main mirror of UTl-Antu. B: 
UTl-Antu main mirror re-aluminisation after the water 
condensation event. C: FORS1 moved from UTl-Antu 
to UT3-Melipal. Plotted zeropoints have been corrected 
for extinction and colour terms using average values (see 
Appendix El . 

In the case of FORS1, observations of photometric 
standard fields (Landolt 1992) are regularly obtained as 
part of the calibration plan; typically one to three fields 
are observed during each service mode night. The pho- 
tometric zeropoints were derived from these observations 
by means of a semi-automatic procedure, assuming con- 
stant extinction coefficients and colour terms. For a more 
detailed discussion on these parameters, the reader is re- 
ferred to Appendix^ where we show that this is a reason- 
able assumption. Figure [3 shows that, with the exception 
of a few cases, Paranal is photometrically stable, being 
the RMS zeropoint fluctuation cr mo =0.03 mag in U and 
er mo =0.02 mag in all other passbands. Three clear jumps 
are visible in Figure[3 all basically due to physical changes 
in the main mirror of the telescope. Besides these sudden 
variations, we have detected a slow decrease in the effi- 
ciency which is clearly visible in the first 10 months and is 
most probably due to aluminium oxidation and dust de- 
position. The efficiency loss appears to be linear in time, 
with a rate steadily decreasing from blue to red passbands, 
being 0.13 mag yr" 1 in U and 0.05 mag yr _1 in /. To al- 
low for a proper compensation of these effects, we have 
divided the whole time range in four different periods, in 
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which we have used a linear least squares fit to the zero 
points obtained in each band during photometric nights 
only. This gives a handy description of the overall sys- 
tem efficiency which is easy to implement in an automatic 
calibration procedure. 

To derive the colour correction included in the cali- 
bration equation one needs to know the sky colours C s k y ■ 
In principle C s k y can be computed from the instrumen- 
tal magnitudes, provided that the data which correspond 
to the two passbands used for the given colour are taken 
closely in time. In fact, the sky brightness is known to have 
quite a strong time evolution even in moonless nights and 
far from twilight (Walker 1988b, Pilachowski et al. 1989, 
Krisciunas 1990, Leinert et al. 1998) and using magni- 
tudes obtained in different conditions would lead to wrong 
colours. On the other hand, very often FORS1 images are 
taken in rather long sequences, which make use of the same 
filter; for this very reason it is quite rare to have close-in- 
time multi-band observations. Due to this fact and to al- 
low for a general and uniform approach, we have decided 
to use constant sky colours for the colour correction. In 
fact, the colour terms are small, and even large errors on 
the colours produce small variations in the corresponding 
colour correction. 

For this purpose we have used color-uncorrected sky 
brightness values obtained in dark time, at airmass z < 1.3 
and at a time distance from the twilights At tW i > 2.5 
hours and estimated the typical value as the average. The 
corresponding colours are shown in Table where they 
are compared with those obtained at other observatories. 
As one can see, B — V and V — R show a small scatter 
between different observatories, while U — B and V — I 
are rather dispersed. In particular, V — I spans almost a 
magnitude, the value reported for Calar Alto being the 
reddest. This is due to the fact that the Calar Alto sky in 
/ appears to be definitely brighter than in all other listed 
sites. 

Now, given the color terms reported in Table lA~2l the 
color corrections 7 x C s k y turn out to be — 0.02±0.02, 
-0.09±0.02, +0.04±0.01, +0.02±0.01 and -0.08±0.02 in 
U, _B, V , R and / respectively. The uncertainties were 
estimated from the dispersion on the computed average 
colours, which is oc —0.3 for all passbands. We empha- 
size that this rather large value is not due to measurement 
errors, but rather to the strong intrinsic variations shown 
by the sky brightness, which we will discuss in detail later 
on. We also warn the reader that the colour corrections 
computed assuming dark sky colours are not necessarily 
correct under other conditions, when the night sky emis- 
sion is strongly influenced by other sources, like Sun and 
moon. At any rate, colour variations of 1 mag would pro- 
duce a change in the calibrated magnitude of ~0.1 mag in 
the worst case. 

Due to the increased depth of the emitting layers, 
the sky becomes inherently brighter for growing zenith 
distances (see for example Garstang 1989, Leinert et al. 
1998). In order to compare and/or combine together sky 
brightness estimates obtained at different airmasses, one 



Table 3. Number of sky brightness measurements ob- 
tained with FORS1 in U, B, V, R, I passbands from April 
1, 2000 to September 30, 2001. 



Filter 


/* (%) 


N t 


N s 


/. (%) 


U 


1.8 


204 


68 


33.3 


B 


11.3 


479 


434 


90.6 


V 


17.3 


845 


673 


79.6 


R 


27.1 


1128 


1055 


93.5 


I 


42.5 


1783 


1653 


92.7 






4439 


3883 


87.5 



needs to take into account this effect. The law we have 
adopted for the airmass compensation and its ability to 
reproduce the observed data are discussed in Appendix El 
(see Ea. IC.3(I . After including this correction in the cali- 
bration equation and neglecting the error on X, we have 
computed the global RMS error on the estimated sky 
brightness in the generic passband as follows: 

*M ahv * yjs*^ y + < + 7 2 o% + C* <j* + (X 1)2 al(2) 

where a ma , oc , o 7 and <jk are the uncertainties on the 
zero point, sky colour, colour term and extinction coeffi- 
cient respectively. Using the proper numbers one can see 
that the typical expected global error is 0.03-^0.04 mag, 
with the measures in V and R slightly more accurate than 
in U, B and /. 

4. ESO-Paranal night sky brightness survey 

The results we will discuss have been obtained between 
April 1 2000 and September 30 2001, corresponding to 
ESO Observing Periods 65, 66 and 67, and include data 
obtained on 174 different nights. During these eighteen 
months, 4439 images taken in the UBVRI passbands and 
processed by the FORS pipeline were analysed and 3883 
of them (~88%) were judged to be suitable for sky bright- 
ness measurements, according to the criteria we have dis- 
cussed in Patat (J2002J) . The numbers for the different pass- 
bands are shown in Table |31 where we have reported the 
total number N t of examined frames, the number N s of 
frames which passed the tests and the percentage of suc- 
cess f s . As expected, this is particularly poor for the U 
filter, where the sky background level is usually very low. 
In fact, we have shown that practically all frames with a 
sky background level lower than 400 electrons are rejected 
( Patat 2002. Sec. 5). Since to reach this level in the U pass- 
band one needs to expose for more than 800 seconds (see 
Table , this explains the large fraction of unacceptable 
frames. We also note that the number of input frames in 
the various filters reflects the effective user's requests. As 
one can see from Table |3 the percentage of filter usage f t 
steadily grows going from blue to red filters, with R and 
I used in almost 70% of the cases, while the U filter is 
extremely rarely used. 

To allow for a thorough analysis of the data, the sky 
brightness measurements have been logged together with 
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Table 2. Typical night sky broad band colours measured at various observatories. 



Site 


Year 


U -B 


B - V 


V-R 


V-I 


Reference 


Cerro Tololo 


1987-8 


-0.7 


+0.9 


+0.9 


+1.9 


Walker (1987, 1988a) 


La Silla 


1978 




+1.1 


+0.9 


+2.3 


Mattila et al. (1996) 


Calar Alto 


1990 


-0.5 


+1.1 


+0.9 


+2.8 


Leinert et al. (1995) 


La Palma 


1994-6 


-0.7 


+0.8 


+0.9 


+1.9 


Benn & Ellison (1998) 


Paranal 


2000-1 


-0.4 


+1.0 


+0.8 


+1.9 


this work 




10 2.Q 3 



Fig. 3. Distribution of telescope pointings in Alt-Az co- 
ordinates. The lower left insert shows the airmass distri- 
bution. 

a large set of parameters, some of which are related to the 
target's position and others to the ambient conditions. 
The first set has been computed using routines adapted 
from those coded by J. Thorstensen 1 and it includes av- 
erage airmass, azimuth, galactic longitude and latitude, 
ecliptic latitude and hclio-ccliptic longitude, target-moon 
angular distance, moon elevation, fractional lunar illumi- 
nation (FLI), target-Sun angular distance, Sun elevation 
and time distance between observation and closest twi- 
light. Additionally we have implemented two routines to 
compute the expected moon brightness and zodiacal light 
contribution at target's position. For the first task we have 
adopted the model by Krisciunas & Schaefer l|1991fl and 
its generalisation to UBVRI passbands (Schaefer 1998), 
while for the zodiacal light we have applied a bi-lincar in- 
terpolation to the data presented by Lcvasseur-Regourd 
& Dumont (1980). The original data is converted from 
Sio{V) to cgs units and the UBRI brightness is com- 
puted from the V values reported by Levasseur-Regourd 
& Dumont assuming a solar spectrum, which is a good 

1 The original routines by J. Thorstensen are available at the 
following ftp server: ftp://iraf.noao. edu/contrib/skycal.tar.Z 



approximation in the wavelength range 0.2— 2/im (Leinert 
et al. 1998). For this purpose we have adopted the Sun 
colours reported by Livingston {2D00), which turn into 
the following U, B, R and / zodiacal light intensities nor- 
malised to the V passband: 0.52, 0.94, 0.77 and 0.50. For 
a derivation of the conversion factor from Sio(V) to cgs 
units, see Appendix 151 

The ambient conditions were retrieved from the VLT 
Astronomical Site Monitor (ASM, Sandrock et al. 2000). 
For our purposes we have included air temperature, rel- 
ative humidity, air pressure, wind speed and wind direc- 
tion, averaging the ASM entries across the exposure time. 
Finally, to allow for further quality selections, for each 
sky brightness entry we have logged the number of sub- 
windows which passed the A-test (see Patat 2002 ) and the 
final number n g of selected sub-windows effectively used 
for the background estimate. 

Throughout this paper the sky brightness is expressed 
in mag arcsec~ 2 , following common astronomical practice. 
However, when one is to correct for other effects (like zo- 
diacal light or scattered moon light), it is more practical 
to use a linear unit. For this purpose, when required, we 
have adopted the cgs system, where the sky brightness 
is expressed in erg s _1 cm~ 2 A -1 sr _1 . In these units 
the typical sky brightness varies in the range 10~ 9 — 10 -6 
(see also Appendix IBl. It is natural to introduce a surface 
brightness unit (sbu) defined as 1 sbu=10 -9 erg s _1 cm~ 2 
A -1 sr _1 . In the rest of the paper we will use this unit to 
express the sky brightness in a linear scale. 

Due to the large number of measurements, the data 
give a good coverage of many relevant parameters. This 
is fundamental, if one is to investigate possible dependen- 
cies. In the next sub-sections we describe the statistical 
properties of our data set with respect to these parame- 
ters, whereas their correlations are discussed in Sees. 15181 

4.1. Coordinate distribution 

The telescope pointings are well distributed in azimuth 
and elevation, as it is shown in Figure[3] In particular, they 
span a good range in airmass, with a few cases reaching 
zenith distances larger than 60° . Due to the Alt-Az mount 
of the VLT, the region close to zenith is not observable, 
while for safety reasons the telescope does not point at 
zenith distances larger than 70°. Apart from these two 
avoidance areas, the Alt-Az space is well sampled, at least 
down to zenith distances Z=50°. At higher airmasses the 
western side of the sky appears to be better sampled, due 
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Fig. 4. Distribution of telescope pointings in galactic (left 
panel) and helio-ecliptic (right panel) coordinates. The 
two histograms show the distribution of galactic and eclip- 
tic latitudes. 



to the fact that the targets are sometimes followed well 
after the meridian, while the observations tend to start 
when they are on average at higher elevation. 

Since sky brightness is expected to depend on the ob- 
served position with respect to the Galaxy and the Ecliptic 
(see Leinert et al. 1997 for an extensive review), it is in- 
teresting to see how our measurements are distributed in 
these two coordinate systems. Due to the kind of scientific 
programmes which are usually carried out with FORS1, 
we expect that most of the observations are performed 
far from the galactic plane. This is confirmed by the left 
panel of Figure 0] where we have plotted the galactic co- 
ordinates distribution of the 3883 pointing included in our 
data set. As one can see, the large majority of the points 
lie at \b\ >10°, and therefore the region close to the galac- 
tic plane is not well enough sampled to allow for a good 
study of the sky brightness behaviour in that area. 

The scenario is different if we consider the helio-ecliptic 
coordinate system (Figure 0] right panel) . The observa- 
tions are well distributed across the ecliptic plane for 
| A — A0I >60° and (3 <+30°, where the contribution 
of the zodiacal light to the global sky brightness can 
be significant. As a matter of fact, the large majority 
of the observations have been carried out in the range 
—30° < (3 < +30°, where the zodiacal light is rather im- 
portant at all helio-ecliptic longitudes. This is clearly vis- 
ible in the upper panel of Figure |SJ where we have over 
imposed the telescope pointings on a contour plot of the 
zodiacal light V brightness, obtained from the data pub- 
lished by Levasseur-Regourd & Dumont (1980). 




so ' nc 

my ( de g) 

Fig. 5. Upper panel: distribution of telescope pointings in 
helio-ecliptic coordinates. Over imposed is a contour plot 
of the zodiacal light V brightness at the indicated levels 
expressed in sbu (1 sbu = 1CP 9 erg s _1 cm~ 2 A -1 sr _1 ). 
Original data are from Levasseur-Regourd & Dumont 
(1980). Lower panel: Zodiacal light V brightness profiles 
at four different ecliptic latitudes expressed in sbu (left 
scale) and mag arcsec -2 (right scale). The brightness in- 
crease seen at | A — A© | >150° is the so-called Gegenschein. 
The horizontal dotted line is placed at a typical V global 
sky brightness during dark time (21.6 mag arcsec -2 ). 



We note that the wavelength dependency of the zodia- 
cal light contribution is significant even within the optical 
range. In particular it reaches its maximum contribution 
in the B passband, where the ratio between zodiacal light 
and typical dark time sky flux is always larger than 30%. 
On the opposite side we have the / passband, where for 
I A — Aq| >80° the contribution is always smaller than 
30% (see also O'Connell 1987). Now, using the data from 
Levasseur-Regourd & Dumont (1980) and the typical dark 
time sky brightness measured on Paranal, we can estimate 
the sky brightness variations one expects on the basis of 
the pure effect of variable zodiacal light contribution. As 
we have already mentioned, the largest variation is ex- 
pected in the B band, where already at |A — A0|=9O° the 
sky becomes inherently brighter by 0.4-0.5 mag as one 
goes from \(3\ >60° to \(3\ =0°. This variation decreases to 
~0.15 mag in the / passband. 

Due to the fact that the bulk of our data has been 
obtained at |/3| <30°, our dark time sky brightness esti- 
mates are expected to be affected by systematic zodiacal 
light effects, which have to be taken into account when 
comparing our results with those obtained at high ecliptic 
latitudes for other astronomical sites (see Sec.EJ. 
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4.2. Moon Contribution 

Another relevant aspect that one has to take into account 
when measuring the night sky brightness is the contribu- 
tion produced by scattered moon light. Due to the sci- 
entific projects FORS1 was designed for, the large ma- 
jority of observations are carried out in dark time, either 
when the fractional lunar illumination (FLI) is small or 
when the moon is below the horizon. Nevertheless, accord- 
ing to the user's requirements, some observations are per- 
formed when moon's contribution to the sky background 
is not negligible. To evaluate the amount of moon light 
contamination at a given position on the sky (which de- 
pends on several parameters, like target and moon eleva- 
tion, angular distance, FLI and extinction coefficient in 
the given passband) we have used the model developed 
by Krisciunas and Schaefer (1991), with the double aim 
of selecting those measurements which are not influenced 
by moonlight and to test the model itself. We have forced 
the lunar contribution to be zero when moon elevation is 
h m < —18° and we have neglected any twilight effects. On 
the one hand this has certainly the effect of overestimat- 
ing the moon contribution when —18° < h m <0°, but on 
the other hand it puts us on the safe side when selecting 
dark time data. 

As expected, a large fraction of the observations were 
obtained practically with no moon: in more than 50% of 
the cases moon's addition is from 10 _1 to 10 -3 the typical 
dark sky brightness. Nevertheless, there is a substantial 
tail of observations where the contamination is relevant 
(200—400 sbu) and a few extreme cases were the moon 
is the dominating source (>600 sbu). This offers us the 
possibility of exploring both regimes. 

4.3. Solar activity 

We conclude the description of the statistical properties 
of our data set by considering the solar activity during the 
relevant time interval. As it has been first pointed out by 
Lord Rayleigh l|1928f) . the airglow emission is correlated 
with the sunspot number. As we have seen in Sec. ^ this 
has been confirmed by a number of studies and is now a 
widely accepted effect (sec Walker 1988b and references 
therein). As a matter of fact, all measurements presented 
in this work were taken very close to the maximum of 
sunspot cycle n. 23, and thus we do not expect to see 
any clear trend. This is shown in Figure El where we have 
plotted the monthly averaged Penticton-Ottawa solar flux 
at 2800 MHz (Covington 1969) 2 . We notice that the so- 
lar flux abruptly changed by a factor ~2 between July 
and September 2001, leading to a second maximum which 
lasted roughly two months at the end of year 2001. This 
might have some effect on our data, which we will discuss 
later on. 



2 The data are available in digital form at the following web 
site: http: //www . drao .nrc . ca/ icarus/www/ archive .html 
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Fig. 6. Penticton-Ottawa Solar flux at 2800 MHz 
(monthly average). The time range covered by the data 
presented in this paper is indicated by the thick line. The 
upper left insert traces the solar flux during the last six 
cycles. 



5. Dark time sky brightness at Paranal 

Due to the fact that our data set collects observations 
performed under a wide range of conditions, in order to 
estimate the zenith sky brightness during dark time it is 
necessary to apply some selection. For this purpose we 
have adopted the following criteria: photometric condi- 
tions, airmass X <1.4, galactic latitude |6| >10°, time 
distance from the closest twilight Attwi >1 hour and no 
moon (FLI=0 or hu < —18°). Unfortunately, as we have 
mentioned in Sec. 14.11 very few observations have been 
carried out at |/3| >45° and hence we could not put a 
very stringent constraint on the ecliptic latitude, contrary 
to what is usually done (see for example Benn & Ellison 
1998). To limit the contribution of the zodiacal light, we 
could only restrict the range of helio-ecliptic longitude 
(| A — A© | >90°). The results one obtains from this se- 
lection are summarized in Tableland Figure where we 
have plotted the estimates of the sky brightness at zenith 
as a function of time. Once one has accounted for the zo- 
diacal light bias (see below) , the values are consistent with 
those reported for other dark sites; in particular, they are 
very similar to those presented by Mattila et al. (1996) 
for La Silla, which were also obtained during a sunspot 
maximum (February 1978). As pointed out by several au- 
thors, the dark time values show quite a strong dispersion, 
which is typically of the order of 0.2 mag RMS. Peak to 
peak variations in the V band are as large as 0.8 mag, 
while this excursion reaches 1.5 mag in the / band. 
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Fig. 7. Zenith corrected sky brightness measured at 
Paranal during dark time (thick dots) from April 1st, 
2000 to September 30th, 2001. The selection criteria are: 
\b\ >10°, |A - A© | >90°, At tm >1 hour, FLP=0 or h m < 
— 18°. Thin dots indicate all observations (corrected to 
zenith). The horizontal dotted lines are positioned at the 
average values of the selected points. The histograms trace 
the distribution of selected measurements (solid line) and 
all measurements (dotted line), while the vertical dashed 
lines are placed at the average sky brightness during dark 
time. 

Table 4. Zenith corrected average sky brightness dur- 
ing dark time at Paranal. Values are expressed in mag 
arcsec -2 . Columns 3 to 7 show the RMS deviation, min- 
imum and maximum brightness, number of data points 
and expected average contribution from the zodiacal light, 
respectively. 
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In Sec. 14. II we have shown that the estimates presented 
in Tableware surely influenced by zodiacal light effects of 
low ecliptic latitudes. To give an idea of the amplitude of 
this bias, in the last column of Table 01 we have reported 
the correction AmzL one would have to apply to the av- 
erage values to compensate for this contribution. This has 
been computed as the average correction derived from the 
data of Levasseur-Regourd & Dumont l|1980f) . assuming- 
typical values for the dark time sky brightness: as one can 
see, AmzL is as large as ~0.3 mag in the B passband. 



Fig. 8. B, V and R dark time sky brightness variations 
as a function of ecliptic latitude. The solid lines trace the 
behaviour expected from Levasseur-Regourd & Dumont 
(1980) data for the different passbands. 

The sky brightness dependency on the ecliptic latitude 
is clearly displayed in Figure|Hl where we have plotted the 
deviations from the average sky brightness (cf. Table 0} 
for B, V and R passbands, after applying the correction 
Attizl- We have excluded the / band because it is heavily 
dominated by airglow variations, which completely mask 
any dependency from the position in the helio-ecliptic co- 
ordinate system; the U data were also not included due 
to the small sample. For comparison, in the same figure 
we have over imposed the behaviour expected on the ba- 
sis of Levasseur-Regourd & Dumont i|1980[l data, which 
have been linearly interpolated to each of the positions 
(A — Xq , f}) in the data set. As one can see, there is a rough 
agreement, the overall spread being quite large. This is 
visible also in a similar plot produced by Benn & Ellison 
(1998, their Figure 10) and it is probably due to the night- 
to-night fluctuations in the airglow contribution. 

6. Sky brightness variations during the night 

In the literature one can find several and discordant results 
about the sky brightness variations as a function of the 
time distance from astronomical twilight. Walker fl988b) 
first pointed out that the sky at zenith gets darker by 
~0.4 mag arcsec -2 during the first six hours after the end 
of twilight. Pilachowski et al. ( 1989) found dramatic short 
time scale variations, while the steady variations were at- 
tributed to airmass effects only (see their Figure 2). This 
explanation looks indeed reasonable, since the observed 
sky brightening is in agreement with the predictions of 
Garstang <|1989[) . Krisciunas (1990, his Figure 6) found 
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that his data obtained in the V passband showed a de- 
crease of ~0.3 mag arcsec -2 in the first six hours after 
the end of twilight, but he also remarked that this effect 
was not clearly seen in B. Due to the high artificial light 
pollution, Lockwood, Floyd & Thompson 1|1990[) tend to 
attribute the nightly sky brightness decline they observe at 
the Lowell Observatory, to progressive reduction of com- 
mercial activity. 

Walker's findings were questioned by Leinert et al. 
( f 995) and Mattila et al. I|1996ll . who state that no indica- 
tions for systematic every-night behaviour of a decreasing 
sky brightness after the end of twilight were shown by their 
observations. Krisciunas l|1997D notes that, on average, the 
zenith sky brightness over Mauna Kea shows a not very 
convincing sky brightness change of 0.03 mag hour -1 . On 
the other hand he also reports cases where the darkening 
rate was as large as 0.24 mag hour -1 and discusses the 
possibility of a reverse Walker effect taking place during 
a few hours before the beginning of morning twilight. 

Leinert et al. (1998) touch this topic in their exten- 
sive review, pointing out that this is an often observed 
effect due to a decreasing release rate of the energy stored 
in the atmospheric layers during day time. Finally, Bonn 
& Ellison ( 199"5|l do not find any signature of steady sky 
brightness variation depending on the time distance from 
twilights at La Palma, and suggest that the effect observed 
by Walker is due to the variable contribution of the zo- 
diacal light, a hypothesis already discussed by Garstang 
( 1997). A further revision of Walker's findings is presented 
here in Appendix [D] where we show that the effect is sig- 
nificantly milder than it was thought and probably influ- 
enced by a small number of well sampled nights. 

We have performed an analogous analysis on our data 
set, using only the measurements obtained during dark 
time and correcting for differential zodiacal light contri- 
bution. Since the time range covered by our observations 
is relatively small with respect to the solar cycle, we do 
not expect the solar activity to play a relevant role, and 
hence we reckon it is reasonable not to normalise the 
measured sky brightness to some reference time. This op- 
eration would be anyway very difficult, due to the vast 
amount of data and the lack of long time series. The re- 
sults are presented in Figure where we have plotted the 
sky brightness vs. time from evening twilight, At etW i, for 
B,V,R and I passbands. Our data do not support the ex- 
ponential drop seen by Walker l|198 8b) during the first 4 
hours and confirm the findings by Leinert et al. (1995), 
Mattila et al. fiGGty and Benn & Ellison lfP9~9"%jl . This is 
particularly true for V and R data, while in B and espe- 
cially in I one might argue that some evidence of a rough 
trend is visible. As a matter of fact, a blind linear least 
squares fit in the range < At e t W i < 6 gives an average 
slope of 0.04±0.01 and 0.03±0.01 mag hour -1 for the two 
passbands respectively. Both values are a factor of two 
smaller than those found by Walker |1988b| but are con- 
sistent, within the quoted errors, with the values we found 
revising his original data (see Appendix lD)l . However, the 
fact that no average steady decline is seen in V and R 
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Fig. 9. Dark time Paranal night sky brightness, corrected 
for zodiacal light contribution, as a function of time dis- 
tance from evening twilight. The vertical dotted lines in- 
dicate the shortest and longest night (7.4 and 10.7 hours 
respectively, astronomical twilight to twilight), while the 
dashed horizontal line is placed at the average value in 
each passband. 

casts some doubt on the statistical significance of the re- 
sults one gets from B and I data. This does not mean 
that on some nights very strong declines can be seen, as 
already pointed out by Krisciunas Ijl997|l . Our data set 
includes several such examples, but probably the most in- 
teresting is the one which is shown in Figure ITU1 where we 
have plotted the data collected on five consecutive nights 
(2000 April 3-7). As one can see, the I data (upper panel) 
show a clear common trend, even though segments with 
different slopes are present and the behaviour shown to- 
wards the end of the night during 2000 April 6 is opposite 
to that of 2000 April 5. This trend becomes less clear in 
the R passband (middle panel) and it is definitely not vis- 
ible in V (lower panel), where the sky brightness remains 
practically constant for about 6 hours. Unfortunately, no 
B data are available during these nights. 

A couple of counter-examples are shown in Figure ITT1 
the upper and lower panels show two well sampled time 
series obtained on the same sky patch, which show that 
during those nights the sky brightness was roughly con- 
stant during the phase where the Walker effect is expected 
to be most efficiently at work. Instead of a steady decline, 
clear and smooth sinusoidal fluctuations with maximum 
amplitudes of ~0.1 mag and time scales of the order of 0.5 
hours are well visible. Finally, to show that even mixed be- 
haviours can take place, in the central panel we have pre- 
sented the R data collected on 23-02-2001, when a number 
of different sky patches was observed. During that night, 
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Fig. 10. Time sequences collected on April 2-7, 2000. The 
data have been corrected for airmass and differential zodi- 
acal light contribution. The vertical dotted line is placed 
at the beginning of morning astronomical twilight. 

the sky brightness had a peak-to-peak fluctuation of ~0.7 
mag and showed a steady increase for at least 4 hours. 
To conclude, we must say that we tend to agree with 

(Ml 



Leinert et al. ( 1995J) that the behaviour shown during sin- 
gle nights covers a wide variety of cases and that there 
is no clear average trend. We also add that mild time- 
dependent effects cannot be ruled out; they are probably 
masked by the much wider night-to-night fluctuations and 
possibly by the patchy nature of the night sky even during 
the same night. 

7. Testing the moon brightness model 

As we have mentioned in Sec. 14. 21 some data have been col- 
lected when the moon contribution to the sky brightness 
is conspicuous and this offers us the possibility of directly 
measuring its effect and comparing it with the model by 
Krisciunas & Schaefer (1991) which, to our knowledge, is 
the only one available in the literature. 

To estimate the fraction of sky brightness generated by 
scattered moon light, we have subtracted to the observed 
fluxes the average values reported in Table 0] for each 
passband. The results are presented in the lower panel 
of Figure ^1 where we have plotted only those data for 
which the observed value was larger than the dark time 
one. As expected, the largest deviations are seen in B, 
where the sky brightness can increase by about 3 mag at 
10 days after new moon, while in /, at roughly the same 
moon age, this deviation just reaches ~1.2 mag. It is in- 
teresting to note that most exposure time calculators for 
modern instruments make use of the function published by 
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Fig. 11. Time sequences collected on 19-12-2000 (I), 23- 
02-2001 (R) and 16-07-2001 (I). The data have been cor- 
rected for airmass and differential zodiacal light contribu- 
tion. The vertical dotted line is placed at the beginning of 
morning astronomical twilight. 

Walker ( 1987) to compute the expected sky brightness as 
a function of moon age. As already noticed by Krisciunas 
(1990), this gives rather optimistic estimates, real data 
being most of the time noticeably brighter. This is clearly 
visible in Figure IT51 where we have overplotted Walker's 
function for the V passband to our data: already at 6 days 
past new moon the observed V data (open squares) show 
maximum deviations of the order of 1 mag. These results 
are fully compatible with those presented by Krisciunas 
(1990) in his Figure 8. 

Another weak point of Walker's function is that it has 
one input parameter only, namely the moon phase, and 
this is clearly not enough to predict with sufficient accu- 
racy the sky brightness. This, in fact, depends on a number 
of parameters, some of which, of course, are known only 
when the time the target is going to be observed is known. 
In this respect, the model by Krisciunas & Schaefer (1991 ) 
is much more promising, since it takes into account all rele- 
vant astronomical circumstances. The model accuracy was 
tested by the authors themselves, who reported RMS de- 
viations as large as 23% in a brightness range which spans 
over 20 times the typical value observed during dark time. 

In the upper panel of Figure IT2*1 we have compared our 
results with the model predictions, including B 7 V, R and 
/ data. We emphasise that we have used average values 
for the extinction coefficients and dark time sky bright- 
ness and this certainly has some impact on the computed 
values. On the other hand, this is the typical configura- 
tion under which the procedure would be implemented in 
an exposure time calculator, and hence it gives a realis- 
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Fig. 12. Lower panel: observed sky brightness variation 
as a function of moon age for B,V,R and /. The solid 
line traces the data published by Walker i|1987|) for the 
V passband while the upper scale shows the fractional lu- 
nar illumination. Upper panel: comparison between the 
observed and predicted moon contribution (Krisciunas & 
Schaefer I1991f) . Plotted are only those data points for 
which the global brightness is larger than the typical dark 
time brightness. 

tic evaluation of the model practical accuracy. Figure El 
shows that, even if deviations as large as 0.4 mag are de- 
tected, the model gives a reasonable reproduction of the 
data in the brightness range covered by our observations. 
This is actually less than half with respect to the one en- 
compassed by the data shown in Figure 3 of Krisciunas & 
Schaefer (199IJ, which reach ~8300 sbu in the V band. 

8. Sky brightness vs. solar activity 

As we have mentioned in Sec. 14.31 during the time cov- 
ered by the data presented here, the solar activity had 
probably reached its maximum. To be more precise, since 
the current solar cycle (n. 23) has a double peak structure 
(see Figure |5J), our measurements cover the descent from 
the first maximum and the abrupt increase to the second 
maximum. Mainly due to the latter transition, the solar 
density flux at 10.7 cm in our data set ranges from 1.2 
MJy to 2.4 MJy, the median value being 1.8 MJy Even 
though this is almost half of the full range expected on 
a typical complete 11 years solar cycle (80—250 MJy), a 
clear variation is seen in the same solar density flux range 
from similar analysis performed by other authors (see for 
example Mattila et al. 1996, their Figure 6). In Figure H~3l 
we show the case of the R passband, where we have plot- 
ted the nightly average sky brightness vs. the solar density 



flux measured during the day immediately preceeding the 
observations. A linear least squares fit to the data (solid 
line) gives a slope of 0.14±0.01 mag arcsec -2 MJy -1 , 
which turns into a variation of 0.24±0.11 mag arcsec" 2 
during a full solar cycle. This value is a factor two smaller 
than what has been reported for B, V (Walker 1988b, 
Krisciunas 1990) and uvgyr (Leinert et al. 1995, Mattila 
et al. 1996) for yearly averages and it is consistent with 
a null variation at the 2 sigma level. Moreover, since the 
correlation factor computed for the data in Figure ITTA is 
only 0.19, we think there is no clear indication for a real 
dependency. 

This impression is confirmed by the fact that a simi- 
lar analysis for the B and V passbands gives an extrapo- 
lated variation of 0.08±0.13 and 0.07±0.11 mag arcsec -2 
respectively. These numbers, which are consistent with 
zero, and the low correlation coefficients (0.08 and 0.11 
respectively) seem to indicate no short-term dependency 
from the 10.7 cm solar flux. Similar values are found for 
the I passband (Am=0.22±0.15 mag arcsec -2 ). These re- 
sults agree with the findings by Leinert et al. (1995) and 
Mattila et al. I|1996f) and the early work of Rosenberg & 
Zimmermann l|1967fl . who have shown that the [OI]5577A 
line intensity correlates with the 2800MHz solar flux much 
more strongly using the monthly averages than the nightly 
averages. For all these reasons, we agree with Mattila et 
al. I)1996|l in saying that no firm prediction on the night 
sky brightness can be made on the basis of the solar flux 
measured during the day preceeding the observations, as 
it was initially suggested by Walker l|1988b|) . A possible 
physical explanation for this effect is that there is some 
inertia in the energy release from the layers ionised by the 
solar UV radiation, such that what counts is the integral 
over some typical time scale rather than the instantaneous 
energy input. 

9. Discussion and conclusions 

Besides being the first systematic campaign of night sky 
brightness measurements at Cerro Paranal, the survey 
we have presented here has many properties that make 
it rather unique. First of all, the fact that it is com- 
pletely automatic ensures that each single frame which 
passes through the quality checks contributes to build a 
continously growing sample. Furthermore, since the data 
are produced by a very large telescope, the measurements 
accuracy is quite high when compared to that generally 
achieved in this kind of study, which most of the time make 
use of small telescopes. Another important fact, related to 
both the large collecting area and the use of a CCD detec- 
tor, is that the usual problem of faint unresolved stars is 
practically absent. In fact, with small telescopes, it is very 
difficult to avoid the inclusion of stars fainter than V=13 
in the beam of the photoelectric photometer (see for ex- 
ample Walker 1988b). The contribution of such stars is 
39.1 Sio(V) (Roach & Oordon H373l Table 2-1) which cor- 
responds to about 13% of the global sky brightness. Now, 
with the standard configuration and a seeing of 1", during 
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Table 5. Dark time zenith night sky brightness measured at various observatories (adapted from Benn & Ellison 
1998). Sio.7em is the Penticton-Ottawa solar density flux at 2800 MHz (Covington 1969). 



Site 


Year 


$10. 7cm 
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B 


v 


R 


/ 


R pforon cp 

_L LL.LL.L LULL 






MJy 




mag 


arcsec 








La Silla 


1978 


1.5 




22.8 


21.7 


20.8 


19.5 


Mattila et al. (1996) 


Kitt Peak 


1987 


0.9 




22.9 


21.9 






Pilachowski et al. (1989) 


Cerro Tololo 


1987-8 


0.9 


22.0 


22.7 


21.8 


20.9 


19.9 


Walker (1987, 1988a) 


Calar Alto 


1990 


2.0 


22.2 


22.6 


21.5 


20.6 


18.7 


Leinert et al. (1995) 


La Palma 


1994-6 


0.8 


22.0 


22.7 


21.9 


21.0 


20.0 


Benn & Ellison (1998) 


Mauna Kea 


1995-6 


0.8 




22.8 


21.9 






Krisciunas (1997) 


Par anal 


2000-1 


1.8 


22.3 


22.6 


21.6 


20.9 


19.7 


this work 




Fig. 13. Lower panel: nightly average sky brightness in 
the R passband vs. solar density flux. The filled circles 
represent data taken at MJD>52114, while the solid line 
indicates a linear least squares fit to all data. Upper panel: 
Penticton-Ottawa solar flux at 2800 MHz during the time 
interval discussed in this paper. The open circles indicate 
the values which correspond to the data presented in the 
lower panel and the dotted line is placed at the median 
value for the solar density flux. 

dark time FORS1 can reach a 5cr peak limiting magnitude 
V ~23.3 in a 60 seconds exposure for unresolved objects. 
As the simulations show (see Patat l2"0"0"2")l . the algorithm 
we have adopted to estimate the sky background is prac- 
tically undisturbed by the presence of such stars, unless 
their number is very large, a case which would be rejected 
anyway by the A-test fPatat [2"002[1 . Now, since the typi- 
cal contribution of stars with V >20 is 3.2 Sio(V) (Roach 
& Gordon fl973|l . we can conclude that the effect of faint 
unresolved stars on our measurements is less than 1%. 

Another distinguishing feature is the time coverage. 
As reported by Benn & Ellison i|1998[l . the large major- 
ity of published sky brightness measurements were carried 



Fig. 14. V, R and / (from top to bottom) dark time 
sky brightness measured at Paranal from April 2000 to 
September 2001. For each passband the average value 
(dashed line) and the ±ler interval (dotted lines) are plot- 
ted. The large solid dots are placed at the seasonal average 
values computed within the time interval indicated by the 
horizontal error bars, while the vertical error bars rep- 
resent the corresponding RMS sky brightness dispersion. 
The labels on the lower side indicate the austral astronom- 
ical season: winter (W), spring (Sp), summer (Su) and fall 
(F). The vertical dashed line indicates December 31, 2000. 



out during a limited number of nights (see their Table 1). 
The only remarkable exception is represented by their own 
work, which made use of 427 CCD images collected on 63 
nights in ten years. Nevertheless, this has to be compared 
with our survey which produced about 3900 measurements 
during the first 18 months of steady operation. This high 
time frequency allows one to carry out a detailed analysis 
of time dependent effects, as we have shown in Sec. |H] and 
to get statistically robust estimates of the typical dark 
time zenith sky brightness. 
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The values we have obtained for Paranal are compared 
to those of other dark astronomical sites in Table 03 The 
first thing one notices is that the values for Cerro Paranal 
are very similar to those reported for La Silla, which were 
also obtained during a maximum of solar activity. They 
are also not very different from those of Calar Alto, ob- 
tained in a similar solar cycle phase, even though Paranal 
and La Silla are clearly darker in R and definitely in /. 
All other sites presented in Table 03 have data which were 
obtained during solar minima and are therefore expected 
to show systematically lower sky brightness values. This 
is indeed the case. For example, the V values measured at 
Paranal are about 0.3 mag brighter then those obtained 
at other sites at minimum solar activity (Kitt Peak, Cerro 
Tololo, La Palma and Mauna Kea). The same behaviour, 
even though somewhat less pronounced, is seen in B and 
/, while it is much less obvious in R. Finally, the U data 
show an inverse trend, in the sense that at those wave- 
lengths the sky appears to be brighter at solar minima. 
Interestingly, a plot similar to that of FigureED a lso gives a 
negative slope, which turns into a variation AU—~ 0.7±0.5 
mag arcsec -2 during a full solar cycle. Due to the rather 
large error and the small number of nights (11), we think 
that no firm conclusion can be drawn about a possible 
systematic effect, but we notice that a similar behaviour 
is found by Leinert et al. (199EJ) for the u passband (see 
their Figure 6). Since the airglow in U is dominated by 
the 2 Herzberg bands A 3 E-X 3 S (Broadfoot & Kendall 
1968), the fact that their intensity seems to decrease with 
an increasing ionising solar flux could probably give some 
information on the physical state of the emitting layers, 
where molecular oxygen is confined. 

At any rate, the BVRI Paranal sky brightness will 
probably decrease in the next 5-6 years, to reach its nat- 
ural minimum around 2007. The expected darkening is 
of the order of 0.4-0.5 mag arcsec -2 (Walker 1988b), but 
the direct measurements will give the exact values for this 
particular site. In the next years this survey will provide 
an unprecedented mapping of the dependency from solar 
activity. So far, in fact, this correlation has been investi- 
gated with sparse data, affected by a rather high spread 
due to the night-to-night variations of the airglow (see for 
instance Figure 4 by Krisciunas 1990), which tend to mask 
any other effect and make any conclusion rather uncertain. 

As already pointed out by several authors, the night 
sky can vary significantly over different time scales, fol- 
lowing physical processes that are not completely under- 
stood. As we have shown in the previous section, even the 
daily variations in the solar ionising radiation are not suf- 
ficient to account for the observed night-to-night fluctua- 
tions. Moreover, the observed scatter in the dark time sky 
brightness (see Sec. EJ is certainly not produced by the 
measurement accuracy and can be as large as 0.25 mag 
(RMS) in the I passband; since the observed distribution 
is practically Gaussian (see Figure^}, this means that the 
/ sky brightness can range over ~1.4 mag, even after re- 
moving the effects of airmass and zodiacal light contribu- 
tion. This unpredictable variation has the unpleasant ef- 



fect of causing maximum signal-to-noise changes of about 
a factor of 2. 

Besides these short time scale fluctuations that we have 
discussed in Sec. HO and the long term variation due to 
the solar cycle, one can reasonably expect some effects 
on intermediate time scales. With this respect we have 
computed the sky brightness values averaged over three 
months intervals, centered on solstices and equinoxes. The 
results for V, R and / are plotted in Figure PHI where we 
have used all the available data obtained at Paranal dur- 
ing dark time, with A tW i >0. This figure shows that there 
is no convincing evidence for any seasonal effect, espe- 
cially in the I passband, where all three-monthly values 
are fully consistent with the global average (thick dashed 
line). The only marginal detection of a deviation from 
the overall trend is that seen in R in correspondence of 
the austral summer of year 2000, when the average sky 
brightness turns out to be ~1.3<r fainter than the global 
average value. Even though a decrease of about 0.1 mag 
is indeed expected in the R passband as a consequence 
of the Nal D flux variation (see Roach & Gordon 1973 
and the discussion below), we are not completely sure this 
is the real cause of the observed effect, both because of 
the low statistical significance and the fact that a simi- 
lar, even though less pronounced drop, is seen at the same 
epoch in the V band, where the Nal D line contribution 
is negligible (see Figure^. 

To illustrate how complex the night sky variations can 
be, we present a sequence of four spectra taken at Paranal 
during a moonless night in Figure [TBI starting more than 
two hours after evening twilight with an airmass rang- 
ing from 1.4 to 2.0. For the sake of simplicity we concen- 
trate on the spectral region 5500-6500A, right at the inter- 
section between V and R passbands, which contains the 
brightest optical emission lines and the so called pseudo- 
continuum (see SecQJ. Due to the increasing airmass, the 
overall sky brightness is expected to grow according to 
Eq. E3 which for V and R gives a variation of about 0.2 
mag. These values are in rough agreement with those one 
gets measuring the continuum variation at 5500A (0.13 
mag) and 6400A (0.18 mag). Interestingly, this is not the 
case for the synthetic V and R magnitudes derived from 
the same spectra, which decrease by 0.32 and 0.51 mag 
respectively, i.e. much more than expected, specially in 
the R band. This already tells us that the continuum and 
the emission lines must behave in a different manner. In 
fact, the flux carried by the [OI]5577A line changes by a 
factor 1.9 from the first to the last spectrum, whereas the 
adjacent continuum grows only by a factor 1.1. For the 
Nal D lines, these two numbers are 1.4 and 1.2, still in- 
dicating a dichotomy between the pseudo-continuum and 
the emission lines. But the most striking behaviour is that 
displayed by the [OI]6300,6364A doublet: the integrated 
flux changes by a factor 5.2 in about two hours and can 
be easily identified as the responsible for the brighten- 
ing observed in the R passband. This is easily visible in 
Figure El where the [OI]6300A component surpasses the 
[OI]5577A in the transition from the first to the second 
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Fig. 15. Evolution of the night sky spectrum on February 
25, 2001 in the wavelength range 5500-6500 A. The origi- 
nal 1800 seconds spectra were obtained with FORS1, using 
the standard resolution collimator and a long slit 1" wide 
(see also the caption of Figure^. In each panel the start- 
ing UT time and airmass X are reported. For presentation 
the four spectra have been normalised to the continuum 
of the first one in the region 5600-5800 A. 



Fig. 16. Night sky spectrum obtained at Paranal on 
February 25, 2002 at 04:53 UT (see Figure EJ. Marked 
are the expected positions for the most common lines pro- 
duced by artificial scattered light (upper ticks) and natural 
atmospheric features (lower ticks). The dotted line traces 
part of the spectrum taken during the same night at 02:39 
UT. 



spectrum and keeps growing in intensity in the subse- 
quent two spectra. The existence of these abrupt changes 
is known since the pioneering work by Barbier (057 ) , who 
has shown that [OI]6300,6364A can undergo strong bright- 
ness enhancements over an hour or two on two active re- 
gions about 20° on either side of the geomagnetic equator, 
which roughly corresponds to tropical sites. With Cerro 
Paranal included in one of these active areas, such events 
are not unexpected. A possible physical explanation for 
this effect is described by Ingham ( 1972), and involves the 
release of charged particles at the conjugate point of the 
ionosphere, which stream along the lines of force of the 
terrestrial magnetic field. We notice that in our example, 
the first spectrum was taken about two hours before local 
midnight, at about one month before the end of austral 
summer. This is in contrast with Ingham's explanation, 
which implies that this phenomenon should take place in 
local winter, since in local summer the conjugate point, 
which for Paranal lies in the northern hemisphere, sees 
the sun later and not before, as it is the case during local 
winter. 

Irrespective of the underlying physical mechanism, the 
[OI]6300,6364A line intensity 3 changed from 255 R to 
1330 R; the fact that the initial value is definitely higher 

3 Line intensities are here expressed in Rayleigh (R). See 
Appendix IbI 



than that expected at these geomagnetic latitudes (<50 
R, Roach & Gordon fl973l Figure 4-12) seems to indicate 
that the line brightening had started before our first obser- 
vation. On the other hand, the intensity of the [OI]5577A 
line in the first spectrum is 220 R, i.e. well in agreement 
with the typical value (250 R, Schubert & Walterschcid 
2000). 

The case of Nal D lines is slightly different, since these 
features follow a strong seasonal variation which makes 
them brighter in winter and fainter in summer, the in- 
tensity range being 30-200 R (Schubert & Walterscheid 
2000). This fluctuation is expected to produce a seasonal 
variation with an amplitude of about 0.1 mag in the R 
passband, while in V the effect is negligible. Actually, the 
minimum intensity of this feature can change from site to 
site, according to the amount of light pollution. In fact, 
most of the radiation produced by low-pressure sodium 
lamps is released through this transition. For example, 
Benn & Ellison ( 1998) report for La Palma an estimated 
artificial contribution to the Sodium D lines of about 70 
R. In our first spectrum, the measured intensity is 73 R, a 
value which, together with the epoch when it was obtained 
(end of summer) and the relatively large airmass (A=1.5), 
indicates a very small contribution from artificial illumina- 
tion. However, a firmer limit can be set analysing a large 
sample of low resolution spectra taken around midsum- 
mer, a task which is beyond the purpose of this paper. 
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To search for other possible signs of light pollution, we 
have examined the wavelength range 3500-5500A of the 
last spectrum presented in Figure IT51 which was obtained 
at a zenith distance of about 60° and at an azimuth of 
313°. A number of Hg and Na lines produced by street 
lamps, which are clearly detected at polluted sites, falls in 
this spectral region. As expected, there is no clear trace 
of such features in the examined spectrum; in particular, 
the strongest among these lines, Hgl 4358A, is definitely 
absent. This appears clearly in Figure El where we have 
plotted the relevant spectral region and the expected po- 
sitions for the brightest Hg and Na lines (Osterbrock & 
Martel 1992). In the same figure we have also marked 
the positions of O2 and OH main features. A compari- 
son with the spectra presented by Broadfoot & Kendall 
(1968) again confirms the absence of the Hgl lines and 
shows that almost all features can be confidently iden- 
tified with natural transitions of molecular oxygen and 
hydroxyl. There are probably two exceptions only, which 
happen to be observed very close to the expected positions 
for Nal 4978, 4983A and Nal 5149, 5163A, lines typically 
produced by high pressure sodium lamps (Benn & Ellison 
1998). They are very weak, with an intensity smaller than 
2 R, and their contribution to the broad band sky bright- 
ness is negligible. Nevertheless, if real, they could indicate 
the possible presence of some artificial component in the 
Nal D lines, which are typically much brighter. This can 
be verified with the analysis of a high resolution spectrum. 
If the contamination is really present, this should show up 
with the broad components which are a clear signature of 
high pressure sodium lamps. The inspection of a low air- 
mass, high resolution (R=43000) and high signal-to-noise 
UVES spectrum of Paranal's night sky (Hanuschik et al. 
2003, in preparation) has shown no traces of neither such 
broad components nor of other Nal and Hgl lines. For this 
purpose, suitable UVES observations at critical directions 
(Antofagasta, Yumbes mining plant) and high airmass pe- 
riodically executed during technical nights, would proba- 
bly allow one to detect much weaker traces of light pol- 
lution than any broad band photometric survey. But, in 
conclusion, there is no indication for any azimuth depen- 
dency in our dark time UBVRI measurements. 

There are finally two interesting features shown in 
Figure ED which deserve a short discussion. The first is 
the presence of Call H&K absorption lines, which are 
clearly visible also in the spectra presented by Broadfoot 
& Kendall (1968) and are the probable result of sunlight 
scattered by interplanetary dust (Ingham 1962). This is 
not surprising, since the spectrum of Figure ITUI was taken 
at /?=— 3?5 and A — A Q =139?8, i.e. in a region were the 
contribution from the zodiacal light is significant (see 
Figure^))). The other interesting aspect concerns the emis- 
sion at about 5200A. This unresolved feature, identified 
as NI, is extremely weak in the spectra of Broadfoot & 
Kendall (jl968), in agreement with its typical intensity (1 
R, Roach & Gordon 1973}. On the contrary, in our first 
spectrum (dotted line in Figure EJ it is very clearly de- 
tected at an intensity of 7.5 R and steadily grows until it 
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Fig. 17. Lower panel: peak signal-to-noise ratio measured 
for the same star on a sequence of 150 seconds I im- 
ages obtained with FORS1 on July 16, 2001. Solid and 
dashed lines trace Eq. [3] for U and / passbands respec- 
tively. Middle panel: seeing measured by the Differential 
Image Motion Monitor (DIMM, Sandrock et al. 2000) at 
5500 A and reported to zenith (empty circles); each point 
represents the average of DIMM data over the exposure 
time of each image. The solid circles indicate the image 
quality (FWHM) directly measured on the images. Upper 
panel: sky background (in ke~) measured on each image. 

reaches 32 R in the last spectrum, becoming the bright- 
est feature in this wavelength range. This line, which is 
actually a blend of several very close NI transitions, is 
commonly seen in the Aurora spectrum with intensities of 
0.1-2 kR (Schubert & Walterscheid 2000) and it is sup- 
posed to originate in a layer at 258 km. The fact that 
its observed growth (by a factor 4.3) follows closely the 
one we have discussed for [OI]6300,6364A, suggests that 
the two regions probably undergo the same micro-auroral 
processes. 

Such abrupt phenomena, which make the sky bright- 
ness variations during a given night rather unpredictable, 
are accompanied by more steady and well behaved varia- 
tions, the most clear of them being the inherent brighten- 
ing one faces going from small to large zenith distances. 
In fact, as we have seen, the sky brightness increases at 
higher airmasses, especially in the red passbands, where 
it can change by 0.4 mag going from zenith to airmass 
X—2. For a given object, as a result of the photon shot 
noise increase, this turns into a degradation of the signal- 
to-noise ratio by a factor 1.6, which could bring it below 
the detection limit. Unfortunately, there are two other ef- 
fects which work in the same direction, i.e. the increase of 
atmospheric extinction and seeing degradation. While the 
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former causes a decrement of the signal, the latter tends 
to dilute a stellar image on a larger number of pixels on 
the detector. Combining Eq. IC.3I the usual atmospheric 
extinction law I = Iq lO" - 4 *^ -1 ) and the law which de- 
scribes the variation of seeing with airmass (s = s X - 6 , 
Roddier 1981) we can try to estimate the overall effect on 
the expected signal-to-noise ratio at the central peak of a 
stellar object. After very simple calculations, one obtains 
the following expression: 



SNR{X) 
SNR 



= x~ 



[(l-f) + fX\~$ io-°-^{x-i) (3) 



where the subscript denotes the zenith (X=l) value. 
Eq.|3is plotted in FigureE|for the two extreme cases, i.e. 
U and I passbands. For comparison, we have overplotted 
real measurements performed on a sequence of I images 
obtained with FORS1 on July 16, 2001. Given the fact 
that the seeing was not constant during the sequence (see 
the central panel), Eq. |3 gives a fair description of the 
observed data, which show however a pretty large scat- 
ter. As one can see, the average SNR ratio decreased by 
about a factor of 2 passing from airmass 1.1 to airmass 
1.6. Such degradation is not negligible, specially when one 
is working with targets close to the detection limit. For 
this reason we think that Eq. |31 could be implemented in 
the exposure time calculators together with the model by 
Krisciunas & Schaefer fTSPTI , to allow for a more accurate 
prediction of the effective outcome from an instrument. 
This can be particularly useful during service mode ob- 
servations, when now-casting of sky conditions at target's 
position is often required. 
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All FORS1 images shown in this paper were obtained dur- 
ing Service Mode runs and their proprietary period has expired. 

Appendix A: Extinction Coefficients and Colour 
Terms 

Extinction coefficients and colour terms were computed 
using the method outlined by Harris, Pirn Fitzgerald & 
Cameron Reed ()1991f) . i.e. via a single-step multi-linear 
least squares fit to the data provided by the observation 
of photometric standard star fields (Landolt 1992). This 
procedure, in fact, allows one to get the photometric so- 
lutions using different stars with suitable colour and air- 
mass ranges, without the need for repeating the observa- 
tions at different airmasses. Due to the FORS1 calibration 
plan implemented during the time range discussed in this 
work, this was a mandatory requirement. To get meaning- 
ful results, the method needs a fair number of stars well 



spread in colour and airmass. While the first requirement 
is almost always fulfilled for the Landolt fields suitable 
for CCD photometry, the second is more cumbersome to 
achieve using FORS1 data. In order to overcome this prob- 
lem we have computed the photometric parameters on a 
bi-monthly basis, for those time ranges where at least 3 
stars were observed at airmass larger than 1.6. Harris et 
al. (1991) recommend to use data spanning at least 1 air- 
mass, but we had to relax this constraint in order to get 
a sufficient number of stars. 

Some tests have shown that the available photomet- 
ric data do not allow the computation of the second or- 
der colour term and extinction coefficient, which are usu- 
ally included in the photoelectric photometry solutions 
(see Harris et al. 1991, their Eqs. 2.9). In fact, the intro- 
duction of such terms in the fitting of our data produces 
random oscillations in the solutions without any signif- 
icant decrease in the variance. Both second order coef- 
ficients are accompanied by large errors and are always 
consistent with zero. This clearly means that the accu- 
racy of our measurements is not sufficient to go beyond 
the first order term. For this reason, we have adopted 
M — m = too + 7 x C fitting function for 

the generic passband, where M and m are respectively 
the catalogue and instrumental magnitudes of the stan- 
dard stars, mo is the zeropoint, 7 is the colour term with 
respect to some colour C, k is the extinction coefficient 
and z is the airmass. To eliminate clearly deviating mea- 
surements, we have computed the photometric solutions 
in two steps. We have first used all data to get a starting 
guess from which we have computed the global RMS de- 
viation a. Then we have rejected all data deviating more 
than 1.5 xct and performed again the least squares fit on 
the selected data. 

The results for the extinction coefficients are presented 
in Table lA~Tl For each time range and passband we have 
reported the best fit value of n, the estimated RMS error 
(both in mag airmass -1 ) and the number of data points 
used for the least squares fit. 

The values of k in the various bands show some minor 
fluctuations. Of course, we cannot exclude night-to-night 
variations, which are clearly observed at other sites (see for 
instance Krisciunas 1990). This would require a dedicated 
monitoring, which is beyond the scope of this work and not 
feasible with the available data. Here we can only conclude 
that, on average, the extinction coefficients do not show 
any significant evolution or clear seasonal effects during 
the 20 months covered by the data discussed in this paper. 
Therefore, given also the purpose of this analysis, we have 
assumed that the extinction coefficients are constant in 
time and equal to the average values reported on the last 
row of Table IA~T1 

The computed colour terms are shown Table E3 We 
note that in those cases where the constraints on the air- 
mass range and the number of data-points were not ful- 
filled, we have performed the best fit keeping k constant 
and equal to the average values given in Table lA~Tl This 
prevented the best fit from giving results with no physical 
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Table A.l. Bi-monthly average extinction coefficients (mag airmass" 1 ) measured at ESO-Paranal with FORS1 during 
photometric nights. For each entry the total number of data points (n) and the number of data points with airmass 
larger than 1.6 (n z ) are reported. The last row of the table shows the weighted means and the corresponding RMS 
errors. 
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88 
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3 


0.09±0.01 
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138 


19 


11/12-00 


0.43±0.01 


78 


9 


0.21±0.01 


86 


16 


0.11±0.01 


89 


13 


0.07±0.01 


83 


14 


0.03±0.01 


85 


10 


01/02-01 


0.45±0.02 


122 


24 


0.22±0.02 


154 


29 


0.10±0.01 


163 


33 


0.07±0.01 


148 


28 


0.04±0.01 


145 


29 


05/06-01 


0.40±0.03 


64 


4 


0.19±0.02 


66 


4 


0.09±0.01 


60 


3 


0.06±0.01 


60 


3 


0.01±0.01 


68 


3 


07/08-01 


0.43±0.02 


117 


18 


0.22±0.01 


139 


17 


0.12±0.01 


143 


20 


0.08±0.01 


124 


17 


0.04±0.01 


68 


17 


09/10-01 


0.44±0.02 


125 


3 


0.24±0.01 


143 


3 


O.lliO.01 


141 


4 


0.08±0.01 


120 


4 


0.06±0.01 


130 


5 




0.43±0.02 






0.22±0.02 






0.11±0.01 






0.07±0.01 






0.05±0.02 







meaning due to the poor airmass coverage. As in the case 
of the extinction coefficients, the color terms show fluc- 
tuations which are within the errors, stronger oscillations 
being seen in the U passband, where the accuracy of the 
photometry is lower. 

On the basis of these results we can conclude that there 
are no strong indications for significant colour term evolu- 
tion and it is therefore reasonable to assume that they are 
constant in time. For our purposes we have adopted the 
weighted mean values shown in the last row of Table lA~2l 

Appendix B: Sky brightness units 

Throughout this paper we have expressed the sky bright- 
ness in mag arcsec" 2 . Since other authors have used differ- 
ent units, we derive here the conversions for the most used 
ones. In the cgs system, the sky brightness is expressed in 
erg s" 1 cm" 2 A" 1 sr" 1 . Now, if m s k y ,x is the sky bright- 
ness in a given passband (in mag arcsec -2 ) and mo, a is 
the photometric system zeropoint of that passband, the 
conversion is obtained as follows: 

B x (cgs) = io-°- 4 (™^>*- m o.A-26.573) (B.l) 

where the constant in the exponent accounts for the 
fact that 1 arcsec 2 corresponds to 2.35xl0~ n sr. For ex- 
ample, if we use the constant mo for the Johnson system 
V band given by Drilling & Landolt J2000), for a typical 
night sky brightness m s k y ,v=21.6 mag arcsec -2 , this gives 
By ~3.7xl0 -7 erg s -1 cm" 2 A -1 sr — 1 . For practical rea- 
sons, we have introduced the surface brightness unit sbu 
(see Sec.0J, which is defined as 1 sbu=10 -9 erg s -1 cm" 2 
A" 1 sr -1 . With this setting, the typical V sky brightness 
is 366 sbu. 

Often B\ is expressed in SI units. The conversion from 
cgs units is simple, since 1 erg s -1 cm" 2 A" 1 sr -1 = 10 W 
m" 2 sr" 1 /xm -1 and hence: 

B\(SI) = lQ- 0A ( m "ky,\-m 0]X -29.073) (B.2) 

so that the typical V sky brightness turns out to be 
3.7xl0" 6 W m" 2 sr" 1 /mi -1 . 



Especially in the past, the sky brightness was expressed 
in Sio, which is defined as the number of 10 th magni- 
tude stars per square degree required to produce a global 
brightness equal to the one observed in the given pass- 
band. Since 1 sr corresponds to 3.282xl0 3 square degrees, 
the conversion is given by the following equation: 

1 Slo(A) = lQ-0.i{l-21-m ,x) erg s -l cm -2 £-1 Sr -1( B .3) 

Again, using the constant given by Drilling & Landolt 
(2000) for the V band, one obtains 1 S W (V) =1.24xl0" 9 
erg s" 1 cm" 2 A" 1 sr" 1 = 1.24 sbu. For U, B, R and I 
the multiplicative constant is 1.38, 2.11, 0.57 and 0.28, 
respectively. Therefore, the sky brightness expressed in 
Sio and in sbu are of the same order of magnitude in 
all optical broad bands. For example, the typical V sky 
brightness is 295 S W {V). 

Another unit used in the past is the nanoLambert (nL) , 
which measures the perceived surface brightness and cor- 
responds to 3.80 Sio{V). For the V passband this gives 
also 1 nL = 4.72 sbu. 

Finally, to quantify the intensity of night sky emission 
lines, the Rayleigh (R) unit is commonly adopted. It is 
defined as 10 6 /47r photons s" 1 cm" 2 sr" 1 and the con- 
version from erg s" 1 cm " 2 sr" 1 is given by the following 
expression: 

1R = 634.4 A erg s" 1 cm" 2 sr" 1 (B.4) 

where A is expressed in A. The intensity of non- 
monochromatic features, like the pseudo-continuum, is 
generally expressed in R A" 1 . 

For a more thorough discussion of the sky brightness 
units the reader is referred to Leinert et al. 1|1998)) and 
Benn & Ellison ifTMSjl . 

Appendix C: Airmass dependency 

The results obtained by Garstang l|1989f) can be used to 
derive an approximate expression for the sky brightness 
dependency on the zenith distance, as already pointed 
out by Krisciunas & Schaefer (1991). If we assume that 
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Table A. 2. Bi-monthly average colour terms for FORS1. The total number of used data points are the same reported 
in Table |A~T1 while the estimated RMS uncertainty on each entry is 0.01. The last row shows the weighted mean values 
and their estimated RMS errors. 



Time Range 


lU-B 


7s -v 


V 

7s -v 


1V-R 


Iv-i 


03/04-00 


0.09 


-0.07 


0.04 


0.02 


-0.05 


05/06-00 


0.08 


-0.08 


0.04 


0.02 


-0.04 


07/08-00 


0.10 


-0.08 


0.03 


0.00 


-0.05 


09/10-00 


0.08 


-0.07 


0.05 


0.04 


-0.04 


11/12-00 


0.04 


-0.09 


0.04 


0.04 


-0.04 


01/02-01 


0.07 


-0.08 


0.03 


0.03 


-0.05 


03/04-01 


0.07 


-0.08 


0.04 


0.03 


-0.04 


05/06-01 


0.03 


-0.09 


0.05 


0.02 


-0.04 


07/08-01 


0.10 


-0.08 


0.03 


0.02 


-0.05 


09/10-01 


0.08 


-0.08 


0.04 


0.04 


-0.03 




0.07±0.02 


-0.08±0.01 


0.04±0.01 


0.03±0.01 


-0.04±0.01 



a fraction / of the total sky brightness is generated by the 
airglow and the remaining (1 — /) fraction is produced out- 
side the atmosphere (hence including zodiacal light, faint 
stars and galaxies), Garstang's Eq. 29 can be rewritten in 
a more general way as follows: 

b(Z)=b [(l-f) + fX)] (C.l) 

where X, the optical pathlength along a line of sight, 
is given by 

X = (1 - 0.96 sin 2 Z)~ 1/2 (C.2) 

Z being the zenith distance. For the hypothesis on 
which this result is based, the reader is referred to 
Garstang's original pubblication and to the analysis pre- 
sented by Roach & Meinel ( 1955). Here we just recall that 
were obtained assuming that the extra- 
terrestrial fraction (1 — /) comes from infinity, while the 
airglow emission is generated in a layer (the so called van 
Rhjin layer) placed at 130 km from the Earth's surface. 

While the intensity of the extra-terrestrial component 
is independent on the zenith distance, this is not true for 
the airglow, due to the variable depth of the emitting layer 
along the line of sight. This is taken into account by the 
term X, which grows towards the horizon. 

In his model Garstang l|1989H considers three differ- 
ent mechanisms for the background light to reach the ob- 
server once it leaves the van Rhjin layer: direct transmis- 
sion, aerosol scattering and Rayleigh scattering. As a mat- 
ter of fact, the first channel dominates on the others (see 
Figure 5 in Garstang 1989) which, in a first approxima- 
tion, can be safely neglected. One is therefore left with 
the propagation of the sky background light across the 
lower atmosphere, which is described by Garstang's Eq. 30 
through the extinction factor EF. Now, if k is the extinc- 
tion coefficient for a given passband (in mag airmass -1 ), 
we can write EF ~ lo -0 - 4 ^" 1 ), so that the expected 
change in the sky brightness as a function of X is given 
by: 

Am = -2.5 log[(l - /) + / X] + k (X - 1) (C.3) 



Table C.l. Zenith corrections (mag) computed according 
to Ea. IC.3l The values have been calculated using Paranal 
average extinction coefficients k (see Tab lA.l|) and /=0.6. 
The zenith distance Z is expressed in degrees, while X is 
in airmasses. 



z 


X 


AU 


AB 


AV 


AR 


AI 


10.0 


1.01 


0.00 


0.01 


0.01 


0.01 


0.01 


20.0 


1.06 


0.01 


0.03 


0.03 


0.03 


0.04 


30.0 


1.15 


0.03 


0.06 


0.08 


0.08 


0.08 


40.0 


1.29 


0.05 


0.11 


0.14 


0.15 


0.16 


50.0 


1.51 


0.07 


0.18 


0.24 


0.26 


0.27 


60.0 


1.89 


0.08 


0.27 


0.37 


0.40 


0.42 



We notice that using f—1 in this equation gives Eq. 2 
of Krisciunas & Schaefer Ijl991(l , which gives slightly larger 
values for Am than Eq. IC.3I 

The expected behaviour is plotted in Figure lC~Tl where 
we present the predicted sky brightness (upper panel) 
and sky colour variations (lower panel) for the different 
passbands. For / we have assumed a typical value of 0.6 
(see Garstang 1989 and Benn & Ellison 1998). As one 
can see the correction becomes relevant ioiZ >30° and 
reaches values as large as ~0.4 mag at Z=60° (see also 
Table IUt|) . These predictions are consistent with the mea- 
surements obtained by Pilachowski et al. JTS89), Mattila 
et al. lfT9~96). Leinert et al. ([1995) and Benn & Ellison 
( 1998). In particular, Pilachowski et al. ( 1989) have clearly 
detected the sky inherent reddening (see their Figure 3), 
which is in fair agreement with the predictions of Eq. IC.3I 
for the (B — V) colour (see also Figure \CA\ lower panel). 

Our data confirm these findings, as it is shown in 
Figure IG2I where we have plotted Paranal dark time 
measurements for B 7 V, R and / passbands. Even though 
the scattering due to night to night variations is quite 
large, the overall trend is consistent with the predictions 
of Ea. IC.3l The U data, which reach airmass X=\.Z6 only, 
do not show a clear sky brightness increase, and this is 
again compatible with the expected value (~0.05 mag). 

The airmass dependency is better visible when consid- 
ering data obtained during single nights, when the over- 
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Fig. C.l. Upper panel: sky brightness variation (in 
mag arcsec -2 ) as a function of zenith distance expected 
from Eq. IC.3I for /=0.6. For each passband the mean 
Paranal extinction coefficients presented in Table lA~Tl were 
adopted. Lower panel: expected colour variation as a func- 
tion of zenith distance. In both plots the upper scale re- 
ports the optical pathlength X , expressed in airmasses. 

all sky brightness is reasonably stable. This is illustrated 
in Figure IC.3I where we present six very good series all 
obtained in the / passband during moonless nights at 
Paranal. For comparison, we have plotted Eq. IC.3I with 
f=0.7, which reproduces fairly well the observed trend 
from airmass 1.2 to 2.0. As one can see, smooth deviations 
with peak values of ~0.1 mag arcsec -1 are detected within 
each single night, while the range spanned by the zenith 
extrapolated values is as large as ~0.7 mag. The fact that 
a comparably large dispersion is seen in the other filters 
(see. Figure lUT]) suggests that night-to-night variations of 
this amplitude are common to all optical passbands. 

We note that on April 7, 2000 and January 1, 2001 the 
observed slope was higher than predicted by Eq. IC.3I This 
implies that the contribution by the airglow to the global 
sky brightness in those nights was probably larger than 
the value we have used here, i.e. 70%. This hypothesis is 
in agreement with the fact that on 07-04-2000 and 01- 
01-2001 the extrapolated zenith value was 0.4—0.5 mag 
brighter than in all other nights shown in Figure IC. 31 

Appendix D: The Walker effect revisited 

In this section we discuss and re-examine the findings pub- 
lished by Walker (1988b). In Figure ITXT1 we have plotted 
the B data of his Table 1, with the exception of July 8, 
1980 and May 4, 1984 (since they cover a very small time 
range with 2 and 3 points only) and the two data points 
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Fig. C.2. B, V, R and / sky brightness during dark time 
as a function of airmass. The solid line traces Eq. IC.3I 
while the dotted lines are placed at ±0.4 mag from this 
law. 

marked by Walker with an asterisk in the sequences of 
May 11, 1983 and Apr 28, 1987. As one can see, there are 
only two dates in which a clear decreasing trend is visi- 
ble, i.e. May 11, 1986 and April 25, 1987. Linear squares 
fitting to the two data sets give slopes of ^0.08 ±0.02 
mag hour -1 both in B and V filters, which turn into a 
sky brightness decrease of 0.48±0.12 mag during the first 
6 hours of the night. Since Walker's data have been col- 
lected across a full sunspot cycle, possible systematic vari- 
ations due to the solar activity have to be removed. This 
can be achieved shifting all time sequences in order to 
have the same sky brightness at some reference time, an 
operation that also has the effect of correcting for night- 
to-night variations in the overall sky brightness. For this 
purpose we have performed, for each date, a linear least 
squares fit to the data and we have interpolated the re- 
sulting straight line at a time distance from the evening 
twilight At etwi —2 hours. We note that a similar proce- 
dure must have been followed by Walker (1988b), since in 
his Figure 2 all measurements refer to the magnitude at 
the end of evening twilight (At etw i=0) . Even though not 
explicitly mentioned in the paper, this implies that some 
extra/interpolation had to be performed, since estimates 
at At e twi=0 are practically never available. 

To avoid meaningless extrapolations we have disre- 
garded the data from July 28, 1981, which do not cover 
a suitable time range. Following Walker, we also have 
not included the data obtained in 1982, since they were 
most likely affected by the Chinchonal volcano eruption. 
Finally, to account for the differential contribution of zodi- 
acal light, we have corrected Walker's measurements using 
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Fig. C.3. Sky brightness sequences on six different moon- 
less nights at Paranal in the / passband, corrected for zo- 
diacal light contribution. The dotted line connects chrono- 
logically the various data points, while the filled dot indi- 
cates the first observation in each time series. Finally, the 
solid line traces Ea. lUlfl with /=0.7. 



the data by Levasseur-Regourd & Dumont (jQSOJ) and as- 
suming that all observations were carried out at the zenith 
of the observing site. The results one obtains following 
this procedure are presented in Figure lD~2l which shows 
a much less convincing evidence for a systematic trend 
than Walker's Figure 2. In that case in fact, a decrement 
of ^0.4 mag is seen during the first four hours, a be- 
haviour which is definitely not visible in our Figure It). 21 
As a matter of fact, a linear least squares fitting gives a 
rate of 0.02±0.01 and 0.03±0.01 mag hour" 1 for B and 
V respectively, which reduce to 0.01±0.01 and 0.02 ±0.01 
mag hour^ 1 if one excludes the two sequences of May 11, 
1986 and April 25, 1987 (empty circles). From these val- 
ues we can deduce a maximum decrement of 0.1 and 0.2 
mag in B and V respectively during the first 6 hours af- 
ter the end of twilight. Therefore, our conclusion is that 
even though there is an indication for a systematic trend 
in Walker's data, the rate is significantly lower than was 
thought. Moreover, due to the limited number of nights 
and the small time coverage of several of the time series, 
the results one gets may depend on the behaviour recorded 
during a few well sampled nights. For example, all data 
at At e twi >5.5 hours were collected on April 28, 1987. 
This could explain why no other author has found the 
same strong effect when using larger data bases (see for 
instance Figure 9 of Benn & Ellison 1998). 
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Fig. D.l. Zenith sky brightness at San Benito Mountain 
in the B passband from Walker l|1988bjl . For clarity, the 
mean sky brightness has been subtracted to each time 
series. 



< 




a 4 
A^lw, (hours) 

Fig. D.2. Variation in zenith sky brightness at San 
Benito Mountain with time after the end of astronomical 
twilight. Data have been corrected for differential zodiacal 
light contribution and all time series were normalised to 
the interpolated brightness at At etW i = 2 hours. Dotted 
and solid lines trace a linear least squares fitting with and 
without the data of May 11, 1986 and April 25, 1987 re - 
spectively (empty circles). Data are from Walker |l988b). 
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